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RADIO AND X-RAY OBSERVATIONS OF SN 2006jd: 
ANOTHER STRONGLY INTERACTING TYPE Hn 

SUPERNOVA 

Poonam Chandra, 1 Roger A. Chevalier, 2 Nikolai Chugai, 3 Claes Fransson, 4 Christopher M. 
Irwin, 2 Alicia M. Soderberg, 5 Sayan Chakraborti 6 & Stefan Immler 7 ' 8 ' 9 

ABSTRACT 

We report four years of radio and X-ray monitoring of the Type Iln supernova 
SN 2006jd at radio wavelengths with the Very Large Array, Giant Metrewave 
Radio Telescope and Expanded Very Large Array; at X-ray wavelengths with 
Chandra, XMM-Newton and Swift-XRT. We assume that the radio and X-ray 
emitting particles are produced by shock interaction with a dense circumstellar 
medium. The radio emission shows an initial rise that can be attributed to free- 
free absorption by cool gas mixed into the nonthermal emitting region; external 
free-free absorption is disfavored because of the shape of the rising light curves 
and the low gas column density inferred along the line of sight to the emission 
region. The X-ray luminosity implies a preshock circumstellar density ~ 10 6 
cm -3 at a radius r ~ 2 x 10 16 cm, but the column density inferred from the 
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photoabsorption of X-rays along the line of sight suggests a significantly lower 
density. The implication may be an asymmetry in the interaction. The X-ray 
spectrum shows Fe line emission at 6.9 keV that is stronger than is expected for 
the conditions in the X-ray emitting gas. We suggest that cool gas mixed into 
the hot gas plays a role in the line emission. Our radio and X-ray data both 
suggest the density profile is flatter than r~ 2 because of the slow evolution of the 
unabsorbed emission. 

Subject headings: Stars: Mass Loss — Supernovae: General — Supernovae: In- 
dividual: SN 2006jd — hydrodynamics — circumstellar matter — radio continuum: 
general 



1. INTRODUCTION 



Type Iln supernovae (SNe Iln) are a subclass of supernovae ( SNe) that show character 
istic narrow emission lines of hydrogen and helium in their spectra (jSchlegellll990t iFilippenko 
19971 ). These SNe are thought to have dense circumstellar media (CSM), and their high Ha 
and bolomet ric luminosities can be explained by th e shock interaction of SN ejecta with the 
dense CSM f lChugailll990uChugai fe Danzigerlll994t ). Observed late excess e s of infrared (1R ) 



radiation are suggestive of dense circumstellar dust flGerardy et al.l 120021 ; iFox et all 1201 it ). 



Generally SNe Iln show significant heterogeneity in terms of their emission line profiles and 
luminosities, and apparent mass loss history. For example, SN 1998S demonstrates a con- 
nection between the Type IIL and the Type Iln class: the light curve and spectra resemble 
those of the bright Type IIL SN 1979C and yet the Ha line shows Type Iln characteristics 
fTFassia et allboOlh . SN 1994W faded within 130 days, whereas SN 1988Z remained bright 
for a decade after the explosion. The Type Iln SN 2005gl has a possible luminous blue 



variable (LBV) star identified as its progenitor ( IGal-Yam fe Leonard 



relatively rare, comprising < 10% of all core-collapse SNe (ISmartt 



2009: 



2009t). SNe Iln are 



Li et al.ll201lt ). 



Dense circumstellar interaction in SNe Iln is expected to produce optical, radio and 
X-ray emission. Radio emission in SNe Iln is expected to be synchrotron emission, initially 
absorb ed mainly by fr ee-free absorption, while X-ray emission is likely to have a thermal 
origin ( IChevalierl Il982t ) . Because of the high density of the CSM, the emission due to the 
CSM interaction is expected to be high. SNe Iln are relatively X-ray luminous, exceeding 



10 



10 



erg s" 



■g., SN 1988Z dFabian fe Terlevichl[l9~96l ISchlegel fe Petrel 120061 ) 



and SN 1986J f lBregman fe Pildid ll992L Il994l : iTemple et al 



20051) . However, SNe Iln have 



been elusive at radio wavelengths. So far only 10 SNe Iln have been detected in the radio. 
van Dyk et al.l (Il996t ) carried out a search for radio emission from 10 SNe Iln and did not 
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detect any. 



SN 2006jd was discovered on 2006 Oct 12, with an apparent magnitude of 17.2, in the 
galaxy UGC 4179 at a redshift of z = 0.0186 JPrasad fc Lill2006h . T he initial spectrum o f 
SN 2006jd showed features of a Type lib SN, similar to SN 1993J falondin et al-lbood ). 
However, based on Keck spectra at late epochs, it was reclassified as a Type Iln SN (Alex 
Filippenko, priv. comm.). On the basis of a spectrum on 2006 O ct 17.51 UT, we assume 
its explosion date to be 2006 Oct 06.5 UT (IBlondin et al.l 120061 ). It is located at R.A. 
= 08/i02m07.43s, Deck = + q0 o 48'31".5 (equinox 2000.0), which is 22" east and 1".3 south 
of the nucleus of UGC 4179. Ilmmler et akl ( 120071 ) detected X-ray emission from SN 2006jd 
on 2007 Nov 16 with the Swift X-ray Telescope (XRT) in a 2.3 ks exposure. The net count 
rate was (6.3 ± 2.0) x 10 -3 cts s _1 . For an adopted plasma model with a temperature of 
10 keV and a Galactic column density of 4.5 x 10 20 cm -2 , they estimated an unabsorbed 
luminosity of (2.5 ± 0.8) x 10 41 erg s _1 , placing it in the category of highly X-ray luminous 
SNe Iln. Radio emission was detected from the SN on 2007 Nov 21.28 UT in the 5 GHz 
band, on 2007 Nov 26.36 UT in the 8.5 GHz band, and on 2007 Nov 26.38 UT in the 22.5 



GHz band flChandra fc Soderberdl2007l ). 



In this paper, we report 4 years of radio and X-ray monitoring of SN 2006jd and its 
implications for the CSM density and the mass loss rate of the progenitor star. In Section 
El we describe the radio and X-ray observations, and their analysis in Section |3j We discuss 
the main results and present conclusions in Section HI 



2. OBSERVATIONS 

2.1. Radio Observations 

Radio observations of SN 2006jd started on 2007 Nov 21.28 UT with the Very Large 
Array @ (VLA). We observed the SN in the 22.5 GHz (K band), 8.5 GHz (X), 5 GHz (C), 
and 1.4 GHz (L) bands at various epochs, along with a 44 GHz band observation at one 
epoch. For C and X bands, the data were taken in the interferometric mode for an average 
of 30 minutes (including calibrator time), whereas for L and K bands the data were collected 
for 1 hr. The total bandwidth used was 50 MHz. 3C48 was used as the flux calibrator and 
J0739+016 was used for the phase calibration. The VLA observed SN 2006jd over the period 
2007 Nov 21.28 to 2009 Aug 7.75. From 2010 Jun onwards, the observations were made using 



1 The Very Large Array is operated by the National Radio Astronomy Observatory, a facility of the 
National Science Foundation operated under cooperative agreement by Associated Universities, Inc. 
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the Expanded VLA (EVLA). For EVLA observations, 3C48 and J0739+016 were also used 
as bandpass calibrators. The data were analyzed using standard AIPS routines. 

We also observed SN 2006jd with the Giant Meterwave Radio Telescope (GMRT) on 
three occasions between 2009 Oct to 2010 Apr in the 1.3 GHz and 0.61 GHz bands. We 
detected the SN in the 1.3 GHz band at both occasions, whereas observations at 0.61 GHz 
resulted in upper limits. For GMRT observations, 3C48 and 3C147 were used as flux calibra- 
tors and J0739+016 as a phase calibrator. All these calibrators were also used for bandpass 
calibration. In the GMRT observations, the observing time on SN 2006jd was approximately 
3 hours in each observation and a bandwidth of 32 MHz was used. AIPS was used to analyze 
the GMRT data as well. 

Details of all the radio observations are listed in Table [TJ We show the radio light 
curves of SN 2006jd at various frequencies in Figure [Hand radio spectra at various epochs 
in Figure [2j 

2.2. X-ray observations 

SN 2006jd was observed with the Swift onboard X-ray Telescope (XRT) at various 
epochs between 2007 Nov and 2011 Mar. In addition, we also observed it once with the 
Chandra and once with the XMM-Newton X-ray observatories. 

We carried out the XMM-Newton observations under our proposal # 55085 starting 2009 
Apr 07 at 4:36:46 UT and continuing until 16:50:06 UT on 2009 Apr 07. The observations 
were carried out with the EPIC-PN and EPIC-MOS cameras in full frame with thin filter 
mode. The exposures for the EPIC-MOS1 and EPIC-MOS2 were 42.667 ks and 42.672 ks, 
respectively, and for the EPIC-PN, the exposure time was 41.032 ks. The XMM-Newton 
on-board Optical Monitor (OM) also observed the SN in UVW1, UVM2 and UVW2 filters 
for 5, 10 and 14.3 ks, respectively. 

We used the EPIC-PN observation to carry out a detailed spectroscopic analysis. The 
XMM-Newton software SAS was used to extract the spectrum. A total of 1963 counts were 
obtained in the 0.2 — 10 keV range, resulting in a net count rate of (6.19 ± 0.17) x 10" 2 cts 
s -1 . The NASA HEASARC software HEAsofi§ was used to do the spectral analysis. 

We also observed SN 2006jd with Chandra under our proposal # 10500688. The ob- 
servations were carried out using ACIS-S with grating NONE in the VFAINT mode. The 



2 http://heasarc. gsfc.nasa.gov/docs/software/lhcasoft/ 
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observations started on 2009 Sept 14 at 00:03:21 UT for a total exposure of 37.24 ks. A total 
of 888 counts were obtained in the 0.2 — 10 keV range with a count rate of (2.38±0.08) x 10~ 2 
cts s _1 . We extracted the spectrum using CIAO software and used HEAsoft to analyse the 
spectrum. 

Swift-XRT observed SN 2006jd on 18 occasions between 2007 Nov and 2011 Mar. All 
the Swift observations were for exposures less than 10 ks. The observations were carried out 
in the photon-counting (PC) mode. To extract the counts from the Swift observations, we 
combined the near simultaneous observations and estimated the count rates using XIMAGE 
and SOSTA of the HEAsoft package. 

Table [2] gives all the X-ray observations used in this paper. Table |3] lists the Swift 
observations which were combined to obtain the count rates. The table gives the best 
estimates of count rates extracted within the 0.2 — 10 keV range. We also obtained the best 
fit positions of SN 2006jd in various Swift observations and estimated the position error in 
arcseconds from that of the actual SN position. 



3. DATA ANALYSIS AND INTERPRETATION 



3.1. Radio data analysis and interpretation 



SN 2006jd remained bright in all the radio bands over the 4 year span of observations, 
except for two early upper limits in L band and occasional upper limits in the K band at later 
epochs. Figure [1] shows the light curve of the SN in L, C, X and K bands. The SN was in the 
optically thick regime in the L band, whereas it reached the optically thin regime in the K 
band. In the X and the C bands, the SN was close to the peak of the light curve, indicating 
a transition from the optically thick to the optically thin phase. The near-simultaneous 
spectra are plotted from day 408 to day 2000 (see Figure [2]). The radio observations of SN 
2006jd show an evolution from a somewhat positive spectral index t o a negative index (see 
Figure [3]). This type of evolution is commonly observed in radio SNe ( jWeiler et al.ll2002l ) and 
is attributed to the transition from optically thick to optically thin radiation. The transition 
can be clearly seen by comparing progressive spectra from the earlier to the late epochs, 
where the peak is shifting from the higher frequencies to lower frequencies with time. 

The primary mechanisms for absorption in the radio emission are free -free absorption 



by the circumstellar materi al (FFA) and synchrotron self-absorption (SSA. IChevalierlll982 



Chevalier fc Franssonl 120031 ) . The relevance of SSA can be determined from the peak radio 
luminosity and the time o f peak, which c a n be related to the velocity of the radio emitting 
region if SSA dominates (j Chevalier 1998 : Chevalier fc Fransson 2003 ). For SN 2006jd, we 
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find a velocity of 2000 — 3000 km s _1 assuming the SSA absorption mechanism, which is 
smaller than expected in a typical SN and also than the velocities deduced from the X-ray 
emission for SN 2006jd (Section l3.3p . We thus conclude that SSA cannot account for sufficient 
absorption and FFA is likely to dominate, although more complicated models (e.g., CSM 
dumpiness or global asymmetry) might relax constraints on the SSA mechanism imposed 
by the expansion velocity. 

We concentrate on FFA models, allowing for the possibility of a relatively flat density 
profile in SN 2006jd. At the high mass loss rates characteristic of SNe Iln, it is pla usible that 
the density profile is no t that of a steady wind. We fit the external FFA model (I Chevalier 



19821 ; IWeiler et al.ll2002[) to the SN 2006id data with a generalized CSM power law density 
dependence p oc r~ s (JFransson et al.lll996l ): 



F(u,t)=K 1 



5 GHz 



t 



tffa = K 2 



1000 day 

2.1 



V 



5 GHz 



exp^-TpFAj 



t 

1000 day 



where a is the optically thin frequency spectral index, which relates to the electron energy 
index 7 (N(E) oc E~"<) as 7 = 1 — 2a. To determine the value of a, we choose the average of 
spectral indices between K and X band and X and C band in Figure [3] after it has become 
optically thin. This value is a = — 1.04±0.05. Thus we fix a = —1.04 (7 = 3.08) in the radio 
absorption models. The parameter 8 is related to the expansion parameter m (R oc t m , where 
R is the shock radius) and the CSM density power law index s (in p oc r~ s ) as 8 = m(l — 2s). 
The assumption that the energy density in the particles and the fields is proportional to the 
postshock energy density leads to (3 = 3m — (3 — a)(ms + 2 — 2m)/2. Here K\ is the radio flux 
density normalization parameter and Ki is the external FFA optical depth normalization 
parameter. Using these parameters, we fit the SN 2006jd radio data to Eq. [TJ We first fit 
the data for standard s = 2 model in Eq. [TJ which gives a very poor fit (Table HJ). The best 
fit 8 here is 8 = 3.46 ± 0.02 implying m = 1.15. Allowing for a flatter density profile, our 
data is best fit with a value m = 1.09 and s = 1.77 (Table HJ). Dashed lines in Figures [U and 
[2] show the external FFA fit to the radio data for SN 2006jd with s = 1.77. The fit to the 
data is not perfect a nd the value m > 1 is a problem because the model for the emission 
requires that m < 1 (j Chevalier! 1 1 9 821 ) . We explored models with m < 1, but found that they 
did not give an acceptable fit to the data. 

A property of the external absorption model is that there is an initial exponential 
(e~ T ) rise in the flux. The data appear to give more of a power law rise, which is also 
observed in other T ype Iln superno va where the rise of the radio flux was followed: SN 1986J 



(IWeiler et al.lll990h and SN 1988Z (Ivan Dvk et al.lll993l : IWilliams et al.ll2002f ). IWeiler et al. 
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( 119901 ) proposed a model in which thermal absorbing gas is mixed into the synchrotron 



emitting gas, so that the flux takes the form 

t Y f 1 ~ exp(-T intFF A; 



1000 day / \ Ti ntFFA 

WFA = K 3 (^) (joo^ ) , (2) 



where again we allow s ^ 2, so that f3 can be expressed as before. IWeiler et al.l ( 119901 ) assume 
that the internal absorbing gas is homologously expanding with density n oc R~ 3 oc t~ 3m . 
Then t mtFFA n 2 R oc t~ 5m , or 5' = —5m. We did not make that assumption in our fits 
because there are other possibilities for the evolution of the density. The internal absorbing 
gas is likely to be in pressure equilibrium with the surrounding hot gas and photoionized by 
the energetic radiation field. If the gas temperature evolves slowly, we have n oc p cx R 2 ~ s t~ 2 
where p is the pressure, so that 5' = — 4 + m(5 — 2s). In our fits, we allow 5' to be a free 
parameter. 

We considered models that allowed for external as well as internal FFA. However, we 
found that external absorption played a negligible role in the best fit case, while the internal 
absorption by itself provided an excellent fit. The fits are shown in Figures Q] and [2] and the 
model parameters are shown in Table HJ A notable feature of the fits is the small magnitude 
of (3, which is associated with the slow evolution in the optically thin regime that can be 
seen in the late light curve evolution (Figure [TJ). Taking s = 2 for the given (3 implies 
m = 1.19, which is not consistent with the hydrodynamic model. However, a density profile 
with s = 1.5 gives m = 0.9, which is a plausible value. In this model, the value of 5' is 
between the values found in the two density evolution models described above. 

It can be seen from equation (J2J) that in the optically thick limit F(y) oc u 2A+a . The 
value of a is determined by the late optically thin evolution. The fact that the optically 
thick spectrum is approximately reproduced (Figure |2]) provides support for the internal 
FAA model. 

The preferred model has an optical depth of unity at 5 GHz at an age of 1000 days. 
Assuming that the absorbing gas is in pressure equilibrium with the X-ray emitting gas and 
that it has a relatively low temperature (10 4 — 10 5 K), it is possible to estimate the mass 
of the gas. We determined the X-ray emitting gas has a temperature T > 20 keV. If we 
take T « 60 keV, corresponding to a shock velocity v s h = 6700 km s" 1 , the density of the 
X-ray emitting gas is 6 x 10 6 cm" 3 at 1000 days (Section 3.3). In pressure equilibrium, the 
density of cool gas is then n c 3.6 x lO 11 ^" 1 cm -3 , where T4 is the temperature in units 
of 10 4 K. The optical depth is r = k v £, where i is the path length through the cool gas and 
k v = 3.6 x 10~ 27 ri 2 T 4 -3 ^ 2 cm -1 is the free-free absorption coefficient. The absorbing mass 
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is M a pa 4irRlrhn c £, where R c is the typical radius of the absorbing gas, fa is the mean 
particle weight, and the optical depth constraint can be used to determine n c £. The result is 
M a k2x 10~ 8 Tf M , showing that a modest amount of cool gas mixed into the emitting 
region can give rise to the needed absorption. The source of the cool gas is likely to radiative 
cooling of dense gas in the shocked region. 

The models presented here assume a constant spectral index for the radiating electrons. 
Physical effects, such as Coulomb losses at low energies and synchr otron losses at high 
energies, can affect the particle spectrum (IFransson fe Bjornssonlll998l ). Estimates of these 
effects indicate that they might be significant, but a complete discussion is beyond the aims 
of this paper. 



3.2. X-ray data analysis 

3.2.1. XMM-Newton 

The X-ray data were best fit with a thermal plasma model at a temperature above 10 
keV, i.e. in a range where the 0.2 — 10 keV spectral shape is not sensitive to temperature. The 
spectrum also indicates a column density of 1.3 x 10 21 cm -2 , which is significantly larger than 
the Galactic absorption column density of 4.5 x 10 20 cm -2 in that direction. To determine 
the significance of the best fit temperature and the column density, we plot 68%, 90% and 
99% confidence contours of column density Nh versus the plasma temperature kT. While 
the column density is a well constrained parameter, the temperature is unconstrained at the 
upper end (see Figure HJ). To better constrain the temperature, we lowered the temperature 
until the \ 2 became significantly worse and the fit was visibly bad. This temperature is 20 
keV. We thus assign 20 keV as the lower limit to the plasma temperature. 

We have a clear detection of a line between 6.5 — 7.0 keV. The line is best fit at energy 
6.76 keV, which, after redshift correction, corresponds to an energy of 6.89 keV. We associate 
this with the 6.9 keV Fe line, which is due to Fe XXVI. There is also an evidence of a possible 
8.1 keV line, which can be attributed to Ni XXVIII. To determine the significance of these 
lines, we plot confidence contours of their energies versus the respective strengths (FigEJ). 
The Fe 6.9 keV line is detected at a high significance level, whereas the 8.1 keV line is 
significant only at the 68% level. 

Now we attempt to fit the whole spectrum. Since the temperature is very high, we 
explore the possiblity of non-thermal emission. The best fit power law photon index of 
T = 1.22 ± 0.09 is so flat that we consider this model unlikely. If we fit the Mekal model 
with a solar abundance, it does not reproduce the Fe 6.9 keV line. However, a Mekal model 
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with 5 times the solar abundance does fit the spectrum well and reproduces the iron line at 
the correct energy but with a narrower width. To look for the possibility of non-equilibrium 
ionization (NEI), we also fit the XSPEC NEI model to our data. The fit is quite good and it 
reproduces the iron line at the correct energy. The best fit ionization time scale in the NEI 
model is r = 6.28 x 10 11 s cm -3 . Since r = nt (n is the number density and t is the age) and 
the XMM spectrum was taken when the SN age was 908 days (or 7.84 x 10 7 s), this implies 
a number density of n = 7.7 x 10 3 cm~ 3 . This seems low and we will discuss the physical 
viability of the NEI model in Section HJ Finally, we fit the thermal bremsstrahlung model 
with one Gaussian line fit as well as two Gaussians. The fit is not improved by adding the 
Gaussian at 8.1 keV. Best fits are given in Tableland our best fit spectra are in Figure [6j 
We chose the bremsstrahlung model along with the Gaussian to be the best fit representation 
to our data. The equivalent width of the Fe line is 1.4 keV and the flux in the iron line is 
about 10% of the total continuum flux in the 0.2 — 10 keV spectrum. 



3.2.2. Chandra- ACIS 

The Chandra spectrum is best fit with a high temperature thermal plasma and a column 
density of 1.6 x 10 21 cm -2 , which agree with our fits to the XMM-Newton spectrum very 
well. Even though the Chandra spectrum is not as detailed as that of the XMM-Newton, we 
clearly detect the 6.9 keV Fe line, though the 8.1 keV line possibly seen in the XMM-Newton 
spectrum is not present. As in the XMM-Newton spectrum, the temperature is not well 
constrained on the upper end, whereas the column density is well constrained. The iron 6.9 
keV line is also well constrained at the 90% confidence level. 

Because of the high temperature, we again look for the possibility of non-thermal emis- 
sion; a fit yields a power law photon index of V — 1.31 ± 0.14, which again is implausibly 
flat. The non-equilibrium ionization model correctly reproduced the Fe 6.9 keV line but 
again indicates a low density (n = 4x 10 3 cm -3 ). The Mekal model with 5 times solar abun- 
dance does reproduce the Fe 6.9 keV line. We fit a bremsstrahlung model with a gaussian 
at 6.9 keV. The details of the Chandra best fit spectrum are in Table [61 Fig. [6] shows the 
XMM-Newton and Chandra spectra of SN 2006jd. 



3.2.3. Swift-XRT 



The exposures in the Swift observations were not large and carrying out an independent 
spectroscopic analysis was not possible with the small number of counts obtained. Thus we 
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converted the count rates to 0.2 — 10 keV fluxes using a thermal plasma model with a 
temperature of 60 keV and column density of 1.3 x 10 21 cm -2 . 

Table [7J gives the unabsorbed fluxes at the various epochs of observation used in this 
paper. We also convert the X-ray fluxes to luminosities using a UGC 4179 distance of 79 
Mpc. Figure [7] shows the 0.2 — 10 keV X-ray light curve for SN 2006jd. The flux evolution 
is slow and is best fit with a power law index of t~°- 2A±0A2 . 



3.3. X-ray Interpretation 

The XMM-Newton and Chandra X-ray spectra can be fit by either an optically thin 
thermal spectrum, which leads to an electron temperature T e > 20 keV, or a nonthermal 
spectrum with photon index T « 1.2. In the nonthermal interpretation, flux oc z/~ ' 2 (as 
a — 1 — r) would be a surprisingly hard spectrum for either synchrotron or inverse Compton 
emission. In the case of synchrotron radiation, one can show that for a plausible emitting 
volume (determined by the supernova shock velocity), the radiating particle would rapidly 
lose energy to synchrotron radiation, tending to steepen the spectrum and making the flat 
spectrum more unlikely. Thus we discard the nonthermal origin of X-rays and proceed with 
the interpretation for the thermal origin of the X-ray spectrum. 

In the thermal interpretation, the temperature T is assumed to be produced by shock 
heating, with the shock velocity v sh = [IQkT / ^^nip)} 1 / 2 = 6700(fcT/60 keV) 1 / 2 km s" 1 , 
where k is Boltzmann's constant, /i is the mean molecular weight, and m p is the proton 
mass. The shock velocity is closer to expectations for the forward shock than the reverse 
shock wave at this age, alt hough X-ray emission from y oung supernovae is typically attributed 



to the reverse shock wave (jChevalier fc Franssonll2003l ). A high velocity reverse shock wave is 



possible if the SN runs into a dense medium which acts as a wall for the expanding supernova 
gas; such a situation is observed in SN 1987A. However, this situation is expected to produce 
rising fluxes at radio and X-ray wavelengths followed by a turnover. The fluxes in SN 2006jd 
are observed to decline slowly from days 400 through 1700, so we consider a forward shock 
origin to be more likely. 

We found that the temperature of the emitting region is higher than can be measured in 
the 0.2 — 10 keV bandpass of the X-ray detectors, so that a spectral luminosity L u is measured 
rather than a total X-ray luminosity. For bremsstrahlung emission we have L u oc r^V/T 1 ^ 2 , 
where V is the emitting volume, so that L u oc f 2m 0-- s )+ 1 _ if the observed L u oc t~ x , then 
x = 0.24 for the present case (Figure [7]). We have x = 2m(s — 1) — 1 or s = 1 + (1 + 
x)/2m = 1 + 0.62/m. Taking m = 0.9 gives s = 1.7. The radio and X-ray observations 
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thus give s ~ 1. 6, which is similar to the value deduced for SN 2010.il fr o m ob servations at 
a earlier epoch ( Chandra et al.ll2012l ). Using results from iFransson et al.l ( 119961 ) for L v as in 
Chandra et al.l ( 120121 ) the mass loss rate M, normalized to R = 10 15 cm, is M_ 3 /v W 2 ~ 5v^ 6 , 
where M_ 3 = M/(10 -3 M yr" 1 ), v w2 is the preshock wind velocity in units of 100 km s" 1 , 
and v± is the average velocity in units of 10 4 km s" 1 at 10 3 days. At 10 3 days, the preshock 
_1 g cm -3 . The corresponding density in the shocked gas is tih ~ 3 x 10 6 



density is 2 x 10 

-3 



-18„ 



cm 



The observed X-ray spectrum is well fit by a single temperature model, and there is no 
sign of emission from a reverse shock region. One possible reason for the lack of such emission 
is that the supernova interaction has evolved to the point where the reverse shock wave has 
moved back to the center and has weakened. In the model described above, the mass swept 
up by the forward shock is M sw = MR s /v w ~ 3M Q . If the ejecta mass is significantly smaller 
than this, the reverse shock region is expected to decline in luminosity. Another possibility 
is that the reverse shock is a cooling shock and the cool shell absorbs the radiation from the 
reverse shock region. In addition, when the reverse shock is radiative, the luminosity from 
the reverse shock region rises as p w , while that from the forward shock region rises as p 2 . In 
SN 1993J, a hot (80 keV) thermal X-ray spectrum was observed in the first few weeks and 
interpreted as emission from the forward shock r e gion, with emission f rom the reverse shock 
absorbed by a cooling shell (jLeising et al. 1994 ; Fransson et al. 1996 ). The conditions for 
radiative cooling at the reverse shock depend on the density profile of the ejecta, and there 
is a transition to nonradiative evolution at late times. At the high densities inferred for SN 
2006jd, that transition can be delayed. 

Here we have assumed that the hot gas is in thermal equilibrium. The time for the 
electro ns and proton s to come into equilibrium is t ep ~ 3 yr for kT = 50 keV and n e = 10 6 



cm 3 ( iSpitzerl Il962l ). At an age of 1000 days, the electrons and protons are thus close 



to equilibrium by Coulomb interactions and may be brought even closer to equilibrium by 
plasma instabilities. If equilibrium has not been obtained, the shock velocity is somewhat 
higher than what we have deduced here. 

A robust result of fitting the X-ray spectrum is the hydrogen column density Nh to the 
source. Of the Nh = 1.3 x 10 21 cm -2 that is observed, 4.5 x 10 20 cm -2 can be attributed 
to the Galaxy, so there is 8.5 x 10 20 cm -2 left for the supernova host galaxy interstellar 
medium and the supernova circumstellar medium. This column density estimate assumes 
that the absorbing medium has solar abundances, and is cool and not highly ionized. In the 
model for the X-ray emission discussed above, the expected column density in the unshocked 
circumstellar medium is 4 x 10 22 cm -2 , a factor of fifty larger than that inferred from the 
observations. A possible reason is that the circumstellar medium is fully ionized by the X-ray 
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radiation. However, the ionization parameter for this case is £ = L/nr 2 = 50, which indicates 
that there may be som e ionization of CNO elements, but not full ionizati on of the gas 



(IVerner fe Ferlandlll996l ). To check on this result, we ran the CLOUDY code (IFerland et al. 



19981 ) for the parameters discussed above. As expected, we found some ionization of CNO 



elements, but little change in the absorption of the X-ray spectrum. 

The temperature of the circumstellar medium is relevant to the absorption of the radio 
emission from the shocked region. In the CLOUDY simulation we found an electron temper- 
ature somewhat less than 10 5 K, which is co nsistent with the ionization parameter. Using 



equation (20) of I Chevalier fe Franssonl (120031 ) and the circumstellar parameters derived from 



3/2 

the X-ray luminosity, the wind optical depth at 5 GHz on day 1000 is r// = 2.4T e5 , where 
T e5 is the electron temperature in units of 10 5 K. The optical depth is thus consistent with 
that indicated by the observations (Figure [2]). The radio emission is expected to be from the 
high temperature X-ray emitting region and we have seen that the absorbing column to the 
X-ray emitting region is relatively small. Assuming that the inner radius is like that found for 
the X-rays and scaling to the lower column density (upper limit) leads to 77/ = 8 x 10" 4 T e5 3 ^ 2 
at 5 GHz on day 1000. The high ionization parameter in this case gives T e5 > 1, so free-free 
absorption along the line of sight to the X-ray emission fails to produce the observed radio 
absorption. This is consistent with our analysis of the radio emission, which showed that 
external free-free absorption does not lead to good model fits to the data. 

Another robust result from the X-ray spectrum is the finding of an Fe emission line 
feature centered at 6.9 keV rest energy (using z = 0.0186), with an equivalent width of 
1.4 keV. The collisional equilibrium model that fits the continuum emission has a weak Fe 
line (0.1 keV equivalent width) at an energy of 6.9 keV. The model line is much weaker 
than that observed. We considered the possibility of ionization non-equilibrium by using 
the nonequilibrium model in XSPEC. The model depends on the postshock temperature, 
which was fit with 60 keV, and the ionization time Tj = nt, which we varied for various 
fits. We found that for = 6 x 10 11 cm -3 s the model provided a good fit to the data, 
with x 2 /dof = 1.03; both the equivalent width and the line energy were well fit. However, 
considering the age of 1000 days, this implies the density to be ~ 7 x 10 3 cm -3 (Section [2. 2[) . 
much less than the density in the hot gas, ~ 3 x lO 6 ^ 1 cm -3 , needed to produce the 
X-ray luminosity. Raising in the nonequilibrium model, we found that by r« = 5 x 10 12 
cm" 3 s, the line was essentially the same as in equilibrium, showing that any density close to 
that needed for the X-ray luminosity fails to produce the Fe line. The nonequilibrium model 
in XSPEC assumes a Sedov blast wave model for the hydrodynamics, which is different from 
the situation studied here, but we do not expect a significant change in the results. 

The equivalent width of Ka line thus exceeds values predicted by ionization equilibrium 
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models. A scenario that might resolve the issue involves the presence of the cool dense shell 
that is formed from cooled gas at the reverse shock wave. The shell is subject to the Rayleigh- 
Taylor (RT) instability so that the forward p ostshock zone turns into a two-phase mixture 
comp osed of rarefied hot and dense cool gas (IChevalier fe Blondinlll995l ; iBlondin fe Ellison 



200ll ). The observed X-ray continuum with a weak Ka line presumably arises from the 



hot component. The two-phase mixture, however, could enhance the production rate of Ka 
photons because of the efficient recombination of hot Fe +26 ions penetrating the cool phase 
with low temperature electrons. This can be illustrated as follows. The total recombination 
rate in the fc-th phase (k is 'cool' or 'hot') is Rk oc PkOikn e ,k, where pk is the probability for 
the fast Fe +26 ion to get into k-th phase, a^ oc is the recombination coefficient, and n e ^ 
is the electron number density. Assuming pk to be proportional to the volume filling factor 
of the corresponding phase (which requires strong mixing), we find that the ratio of recom- 
bination rates of fast Fe +26 ions in the cool and hot phases is R c h = (M c /M h )(T h /T Cte ff) v , 
where Mk is the mass of the corresponding phase and T Cj6 // = max[T c , Th(m e /56m p )} is the 
effective temperature that determines the relative velocity of cool electrons and hot Fe +26 
ions. Adopting T h = 50 keV, T c = 0.001 keV and r] = 0.7 we obtain R Cjh ^ 2 x 10 3 (M C /M h ), 
which suggests that even a low-mass cool phase M c ~ 10~ 3 M/j could dominate the Ka line 
emission. Keeping in mind that the cool phase does not contribute to the X-ray continuum 
at 6.9 keV, we thus conclude that the proposed mechanism could significantly increase the 
equivalent width of the Ka line compared to a single phase hot plasma. However, it is not 
clear whether the required penetration of Fe ions into the cold gas could take place. A 
magnetic field could inhibit the motion of Fe ions. Also, if magnetic field effects are small, 
there is likely to be an intermediate temperature zone between the hot and cool phases, and 
the Fe +26 ions could recombine there. Our main point is that the cool gas presents some 
possibilities for Fe line emission. We argued in Section 3.1 that the absorption of the radio 
emission may also be related cool gas that is mixed into the shocked, hot gas. 



4. DISCUSSION AND CONCLUSIONS 

The basic interpretation of the X-ray and radio data in terms of standard spherical 
models led to inconsistencies. One is that the column density of matter to the X-ray emission 
is about a factor 50 smaller than that needed to produce the X-ray luminosity. Since the 
absorption column is oc p w and the X-ray luminosity to p^, one possibility is that the 
luminosity is enhanced by dumpiness of the gas. However, dumpiness would raise the shock 
velocity needed to produce a particular X-ray temperature and the temperature is already 
close to that expected behind the forward shock, so this solution does not seem feasible. 
The situation, however, may be different in the case that the circumstellar clouds interact 
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primarily with the cool dense shell. Then the cloud shock velocity is large enough to account 
for the high temperature derived from the observed X-ray spectrum. The straightforward 
explanation is then that there is a global asymmetry in the distribution of the circumstellar 
gas that allows a low column density in one direction, while dense intera ction is taking place 
over m uch of the rest of the solid angle (as viewed from the supernova). IChugai fc Danziger 
( 119941 ) had suggested a scenario with equatorial mass loss to explain the presence of fast 
and intermediate velocity shock fronts in the Type Iln SN 1988Z. Polarization observations 
of Type I In superno vae have sh own evidence for lar ge scale asymme try, e.g., SN 1997 eg 
rtrloffman et al.ll2008h . SN 1998S Leonard et al.lhooof ). and SN 2010jl ( TPatat et al-lboilh . 



In this view, the column density to the radio emission would be small, as for the X-ray 
emission, because they are both from the same region. The external free-free absorption 
model for the radio absorption thus fails, and an absorption mechanism internal to the 
emission is indicated. We have argued that thermal absorption internal to the emitting region 
i s a plausible mech anism, as has previously been proposed for other Type Iln supernovae 
rtWeiler et al1ll99ol ). 



There is evidence that the X-ray properties we have found for SN 2006jd are not unusual 
for a Typ e Iln event. Chandra observations of SN 2 001em and SN 2005dk have shown hard 
emission ( IPooley Sz Lewinl l2004t iPooley et all 120071 ). implying a high temperature. In the 
case of SN 200 lem, the co lumn density is lower than expected from the X-ray luminosity 
( IChugai fc Chevalier! 120061 ). as in SN 2006jd. Another point of comparison is the Type Iln 
SN 2010jl, which ha d an unabsorbed 0.2 — 10 keV luminosity of ~ 7 x 10 41 erg s _1 at ages of 
2 months and 1 yr ( Chandra et alifeoid ). Although the X-ray luminosity is just a factor of 
2 higher than that of SN 2006jd, the absorbing column, 10 24 cm -2 at 2 months and 3 x 10 23 
cm -2 at 1 yr, is much higher. In this case, the observed column is roughly consistent with 
that needed to produce the X-ray luminosity. The implication is that the observed line of 
sight passes through the dense preshock gas for SN 2010jl, but not for SN 2006jd. 

In this paper we have treated the radio and X-ray properties of SN 2006jd. A more com- 
prehensive study would include observations at other wavele ngths. SN 2006jd was detected 



as an infrared source at an age ~ 1150 days (IFox et al.ll201ll ). Its luminosity was ~ 10 erg 



s _1 which is larger than the X-ray luminosity observed at this time. If the interpretation o f 
this emission as from circumstellar dust heated by shock emission is correct (IFox et al .1120 111 ), 
it appears that some of the X-ray emission from the shocked region is shielde d from the ob 



serve r. An infrared/X-ray ratio > 1 has also been observed in other SNe Iln (IGerardy et al. 
2002h . 



Our observations at both radio and X-ray wavelengths imply that the circumstellar 
density profile (p oc r -2 ) is flatter than the s = 2 case that would be expected for a steady 
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wind. We find s = (1.5 — 1.6); deviation from the steady case is plausible for a Type Iln 
supernova because the mass loss may be due to an eruptive event. In ad dition, the optica l 



luminosity evolution of SN 2006jd shows a second peak at ~ 500 days ( jStritzingerl |2012[ ) 
that is indicative of a nonstandard density distribution. 
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Table 1. Radio observations of SN 2006jd 



Tel. 


Date of 


Days since 


Freq. 


Flux density 


rms 




Observation 


explosion 


GHz 


ujy 


ujy 


VLA 


2007 Nov 21.28 


404.74 


8.46 


238 


38 


VLA 


2007 Nov 26.38 


409.84 


4.86 


170 


38 


VLA 


2007 Nov 26.38 


409.84 


22.46 


324 


54 


VLA 


2008 Mar 07.23 


511.69 


8.46 


305 


46 


VLA 


2008 Mar 07.24 


511.70 


4.86 


232 


37 


VLA 


2008 Mar 07.26 


511.72 


22.46 


428 


41 


VLA 


2008 Mar 07.30 


511.76 


1.43 


< 222 


74 


VLA 


2008 May 02.88 


568.34 


43.31 


< 1680 


560 


VLA 


2008 May 08.15 


573.61 


8.46 


525 


50 


VLA 


2008 May 08.17 


573.63 


4.86 


480 


54 


VLA 


2008 May 10.18 


575.64 


1.43 


< 291 


97 


VLA 


2008 May 16.00 


581.46 


8.46 


509 


12 


VLA 


2008 May 16.94 


582.40 


1.43 


209 


63 


VLA 


2008 May 25.96 


591.42 


22.46 


683 


229 


VLA 


2008 Dec 15.38 


794.84 


8.46 


1024 


12 


VLA 


2008 Dec 15.40 


794.86 


4.86 


1035 


51 


VLA 


2008 Dec 15.41 


794.87 


1.43 


250 


65 


VLA 


2008 Dec 15.45 


794.91 


22.46 


624 


90 


VLA 


2009 Feb 03.16 


844.62 


8.46 


1033 


23 


VLA 


2009 Feb 03.22 


844.68 


4.86 


1192 


26 


VLA 


2009 Feb 03.28 


844.74 


1.43 


247 


13 


VLA 


2009 Feb 03.34 


844.80 


22.46 


665 


13 


VLA 


2009 Apr 02.02 


902.48 


8.46 


1069 


71 


VLA 


2009 Apr 02.04 


902.50 


4.86 


1178 


55 


VLA 


2009 Apr 02.06 


902.52 


1.43 


552 


17 


VLA 


2009 Apr 02.10 


902.56 


22.46 


390 


55 


VLA 


2009 Aug 07.67 


1030.13 


8.46 


770 


50 


VT A 
V 1 . . \ 


9000 Ann- 07 fiQ 

zuuy -rvug ui.oy 


1 0*^0 1 ^ 




lOUo 


KO 
DZ 


VLA 


2009 Aug 07.71 


1030.17 


1.43 


743 


123 


VLA 


2009 Aug 07.75 


1030.21 


22.46 


< 645 


215 


GMRT 


2009 Oct 21.00 


1104.46 


1.28 


730 


124 


GMRT 


2010 Feb 09.65 


1216.11 


1.28 


674 


79 


GMRT 


2010 Apr 24.52 


1289.98 


0.61 


< 441 


117 


EVLA 


2010 Jun 06.91 


1333.28 


4.96 


2216 


54 


EVLA 


2010 Jul 25.67 


1382.13 


8.46 


918 


37 


EVLA 


2011 Jul 09.72 


1731.18 


4.50 


1659 


13 


EVLA 


2011 Jul 09.72 


1731.18 


7.92 


1045 


33 


EVLA 


2011 Jul 30.65 


1752.11 


22.46 


< 405 


135 


EVLA 


2011 Jul 31.72 


1753.18 


1.39 


1642 


96 


EVLA 


2012 Apr 01.96 


1998.42 


4.50 


1385 


350 


EVLA 


2012 Apr 01.96 


1998.42 


7.90 


966 


25 


EVLA 


2012 Apr 03.02 


1999.48 


1.39 


1564 


208 


EVLA 


2012 Apr 07.99 


2004.45 


22.46 


468 


66 



Note. - 



Here the upper limits arc 3— a upper limits. 
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Table 2. Details of X-ray observations of SN 2006jd 



Date of 


1\/Ti ccinn 


Inst 


Obs. 


T+YT~inQ1 lT*fi 
UAUU3 LiJ. L- 


Obsci'Vditioii 






ID 


ks 


2007 Nov 16.73 


Swift 


XRT 


00031016001 


2.11 


2007 Nov 21.49 


Swift 


XRT 


00031016002 


2.92 


2007 Nov 22.03 


Swift 


XRT 


00031016003 


1.61 


2007 Nov 22.84 


Swift 


XRT 


00031016004 


3.06 


2007 Dec 18.01 


Swift 


XRT 


00031016006 


9.55 


2008 Jan 15.32 


Swift 


XRT 


00031016007 


9.26 


2008 Feb 20.18 


Swift 


XRT 


00031016008 


3.93 


2008 Feb 21.32 


Swift 


XRT 


00031016009 


6.06 


2008 Sep 09.02 


Swift 


XRT 


00031016010 


1.76 


2008 Sep 10.82 


Swift 


XRT 


00031016011 


2.49 


2008 Sep 11.02 


Swift 


XRT 


00031016012 


7.36 


2009 Apr 07.20 


XMM 


EPIC-PN 


0550850101 


42.92 


2009 Sep 10.12 


Swift 


XRT 


00031016013 


9.79 


2009 Sep 14.00 


Chandra 


ACIS-S 


10127 


37.24 


2011 Mar 08.32 


Swift 


XRT 


00037767006 


1.06 


2011 Mar 08.38 


Swift 


XRT 


00037767004 


0.99 


2011 Mar 08.46 


Swift 


XRT 


00037767005 


1.01 


2011 Mar 09.32 


Swift 


XRT 


00037767007 


3.72 


2011 Mar 10.81 


Swift 


XRT 


00037767008 


3.12 


2011 Mar 12.41 


Swift 


XRT 


00037767009 


2.35 



Table 3. Count rates obtained from Swift-XHT Observations 



ObsID 


Exposure 




RA 


Decl 




Error 


Count Rate 








s 


J2000 


J2000 


tt 






cts/s 








0031016001-004 


9695.659 


OS 


02:07.2 


+00 


18 


31.9 


3.47 


(6.816 


± 


1.100) 


X 




-03 


0031016006 


9547.694 


08 


02:07.1 


+00 


18 


34.0 


5.55 


(6.464 


± 


1.100) 


X 


io- 


-03 


0031016007 


9261.903 


08 


02:07.3 


+00 


18 


33.7 


2.94 


(7.421 


± 


1.200) 


X 


10" 


-03 


0031016008-009 


9991.502 


08 


02:07.5 


+00 


18 


33.3 


2.00 


(6.011 


± 


1.000) 


X 


10" 


-03 


00310160010-012 


11621.214 


08 


02:07.1 


+00 


18 


28.9 


5.59 


(8.094 


± 


1.100) 


X 


10" 


-03 


0031016013 


9750.140 


08 


02:07.4 


+00 


18 


29.5 


2.05 


(4.549 


± 


0.910) 


X 


10" 


-03 


0037767004-009 


12151.953 


08 


02:07.5 


+00 


18 


30.4 


1.52 


(5.608 


± 


0.910) 


X 


io- 


-03 



Table 4. Model fits to SN 2006jd radio data (assuming a = -1.04) 



Model 


x 2 /" 




Par am- 1 








Param-2 


Param-3 




Param-4 


Param-5 


External FFA (s = 2) 


2.41 


K ± 


= (2.27 ± 0.03) 


X 


10 3 


K 2 


= 0.46 ± 0.01 


5 = 3.46 ± 0.02 








External FFA (flatter s) 


2.11 


Kl 


= (2.05 ± 0.03) 


X 




K 2 


= 0.40 ± 0.01 


s = 1.77 ± 0.02 


m 


= 1.09 ± 0.01 




Internal FFA 


1.26 


Kl 


= (2.30 ± 0.04) 


X 


10 3 


K S 


= 1.20 ± 0.04 


= 0.46 ± 0.04 


6' 


= 3.15 ± 0.07 




Internal+External FFA 


1.30 


Kl 


= (2.31 ± 0.04) 


X 


10 3 


K 2 = 


= 0.002 ± 0.003 


K 3 = 1.22 ± 0.05 


6 


= 3.58 ± 0.04 


5' = 3.16 ± 0.07 



Note. — see 
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Table 5. Spectral model fits 


to the XMM-Newton spectrum 




Model 


X 2 /v N H Param-1 


Param-2 Abs. flux Unabs. flux 


Total flux 



PowerLaw 


1.03(123) 


1 


41+0-37 
-0.30 


X 


10 21 


r = 


. ,,,+0.15 
1J1 _0.11 






2.90 


X 


10" 


13 


3.33 


X 


10" 


13 


3.78 X 


10" 


13 


-f- Gaussian 










E - 


6 76+ ' 10 
°- ,D -0.10 


EW = 1. 02+5J 


10 
16 


4.57 


X 


10" 


14 


4.58 


X 


10" 


14 








NEI 


1.09(125) 


1 


99 +0.24 
Zi! -0.18 


X 


10 21 


kT = 


79.65j: 48u2 


r = 6.28i 2 ;«fx 


10 11 


3.25 


X 


10" 


13 


3.59 


X 


10" 


13 


3.59 X 


10" 


13 


Mekal 


1.12(126) 


1 


oc + 0.22 
JO -0.21 


X 


10 21 


kT = 


36.77_ 20 18 






3.09 


X 


10" 


13 


3.48 


X 


10" 


13 


3.48 X 


Ki- 


13 


Bremsstrahlung 


1.02(123) 


1 


9^+0.21 
•^°-Q.23 


X 


10 21 


kT = 


56.76tf 8 ;l 3 






2.91 


X 


10" 


13 


3.26 


X 


10" 


13 


3.64 x 


ll)" 


13 


+ Gaussian 












E = 




EW = 1.38t° 


1 8 
13 


3.75 


X 


10" 


14 


3.76 


X 


10" 


14 









Note. — Here iVff is in cm -- , E and EW are in keV, T in cm - s, and the fluxes are in erg cm - s — . The fluxes mentioned are in the 0.2 — 10.0 kcV 
energy range. The abs. and unabs. fluxes are of that particular component in the model, whereas the total flux is for the total unabsorbed flux of the 
model. 



Table 6. Spectral model fits to the Chandra spectrum 



Unabs. flux 



PowerLaw 


0.90(48) 


, s ,+0.61 
1 ' S ' J -0.53 


X 


10 21 


r 


- 1 34+ - 17 

— J-'^-0 .16 




2.82 X 


10" 


13 


3.24 x 


10~ 


13 


3.44 x 


10~ 


13 


+ Gaussian 








E 


= 6 - 67 -0.21 


EW = 0.73+ 1 ;. 12 


1.98 X 
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Note. — Here is in cm , E and EW are in keV, r in cm s, and the fluxes arc in erg cm s . The fluxes mentioned arc in the 0.2—10.0 keV 

energy range. The abs. and unabs. fluxes arc for that particular component in the model, whereas the total flux is for the total unabsorbed flux of the 



Table 7. Unabsorbed 0.2 — 10 keV X-ray fluxes of SN 2006jd at various epochs 



Date of 


Days since 


Instrument 


Flux 


Luminosity 


Observation 


Explosion 




10 _13 erg/cm 2 /s 


10 41 erg/s 


2007 Nov 16.73-22.84 


403.2 ± 3.1 


Swift-XRT 


4.51 ± 0.73 


3.37 ±0.54 


2007 Dec 18.01 


431.5 


Swift-XRT 


4.28 ±0.73 


3.19 ±0.54 


2008 Jan 15.32 


459.8 


Swift-XRT 


4.91 ± 0.79 


3.67 ±0.59 


2008 Feb 20.18-21.32 


496.2 ± 0.6 


Swift-XRT 


3.98 ± 0.66 


2.97 ±0.49 


2008 Sep 09.02-11.02 


698.5 ± 1.0 


Swift-XRT 


5.35 ±0.73 


4.00 ± 0.54 


2009 Apr 07.20 


907.7 


XMM EPIC-PN 


q f?q+0.33 
°- uo -0.30 


o 70 + 0.25 
z - ' u -0.23 


2009 Sep 10.12 


1063.6 


Swift-XRT 


3.01 ±0.60 


2.25 ± 0.45 


2009 Sep 14.00 


1067.5 


Chandra ACIS-S 


o no + 0.89 

°' oo -0.93 


9 co+0.70 
Z - OZ _0.66 


2011 Mar 08.32-12.41 


1609.8 ± 2.0 


Swift-XRT 


3.71 ±0.60 


2.77 ±0.45 
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Fig. 1 — Radio light curves of SN 2006 in the 1.4 GHz (L), 5 GHz (C), 8.5 GHz (X) and 22.5 
GHz (K) bands. The radio emission in the L band is still rising and the K band emission 
appears to be in the optically thin phase at late times. The radio emission is making a 
transition from optically thick to thin in the X and C bands. Here dashed lines show the 
best fit external absorption model for s = 1.77. Black solid lines are the best fit internal 
absorption model. 
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Fig. 2. — Radio spectra of SN 2006jd at 9 epochs. Even though the spectral evolution is 
slow, a transition from optically thick to an optically thin is clearly seen in the progression 
of the spectra from early to late epochs. Dashed lines show the best fit external absorption 
model for s = 1.77. Black solid lines are the best fit internal absorption model. 
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Fig. 3. — Spectral index (a) evolution for SN 2006jd between various radio frequencies. The 
evolution is from positive index towards negative, typical in core-collapse SNe. 




N H (10 ZZ cm- 2 ) 

Fig. 4. — Confidence contours of column density versus the X-ray emitting plasma tempera- 
ture. The column density is well constrained whereas the temperature is not constrained at 
the higher end. 
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Fig. 5. — Confidence contours of the Fe 6.9 keV line (left panel) and the 8.1 keV line (right 
panel). The iron line is detected at a high significance level whereas detection of the 8.1 keV 
line is not significant above 1-sigma level. 




Fig. 6. — The left panel shows the XMM-Newton spectrum of SN 2006jd and the right panel 
shows the Chandra spectrum of SN 2006jd. Note that the 6.9 keV Fe line is detected in 
both spectra and has 10% of the total continuum flux (estimated from the XMM-Newton 
spectra). 
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Fig. 7. — X-ray light curve of SN 2006jd. The light curve is best fit with a power law of 
index -0.24 ±0.12. 



